A sweep-based low-rank method for the discrete ordinate transport equation

نویسندگان

چکیده

The dynamical low-rank (DLR) approximation is an efficient technique to approximate the solution matrix differential equations. Recently, DLR method was applied radiation transport calculations reduce memory requirements and computational costs. This work extends scheme for time-dependent equation in 2-D 3-D Cartesian geometries with discrete ordinates discretization angle (SN method). reduced system that evolves on a manifold constructed via “unconventional” basis update & Galerkin integrator avoid substep backward time, which could be unstable dissipative problems. resulting preserves information angular direction by applying separate decompositions each octant where intensity has same sign as cosines. Then, sweeps source iteration can efficiently solve this low-rank-SN system. numerical results demonstrate requires less time than solving full rank equations using without losing accuracy.

برای دانلود باید عضویت طلایی داشته باشید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

A uniformly second order numerical method for the one-dimensional discrete-ordinate transport equation and its diffusion limit with interface

Abstract. In this paper, we study a uniformly second order numerical method for the discreteordinate transport equation in the slab geometry in the diffusive regimes with interfaces. At the interfaces, the scattering coefficients have discontinuities, so suitable interface conditions are needed to define the unique solution. We first approximate the scattering coefficients by piecewise constant...

متن کامل

An adaptive Runge-Kutta-Fehlberg method for time-dependent discrete ordinate transport

To my beautiful wife, Ashley. For your unconditional love, faithful encouragement, and resolute support. You have always believed in and inspired the pursuit of my dreams. iv ACKNOWLEDGEMENTS I would first like to express my immense gratitude to my advisor, Dr. Glenn Sjoden. What started as an undergraduate summer research project transformed into my path in graduate school. As both an advisor ...

متن کامل

Discrete-Ordinate Discontinuous Galerkin Methods for Solving the Radiative Transfer Equation

Abstract. The radiative transfer equation (RTE) occurs in a wide variety of applications. In this paper, we study discrete-ordinate discontinuous Galerkin methods for solving the RTE. The numerical methods are formed in two steps. In the first step, the discrete ordinate technique is applied to discretize the integral operator for the angular variable, resulting in a semi-discrete hyperbolic sy...

متن کامل

A Riemannian rank-adaptive method for low-rank optimization

This paper presents an algorithm that solves optimization problems on a matrix manifold M ⊆ Rm×n with an additional rank inequality constraint. The algorithm resorts to well-known Riemannian optimization schemes on fixed-rank manifolds, combined with new mechanisms to increase or decrease the rank. The convergence of the algorithm is analyzed and a weighted low-rank approximation problem is use...

متن کامل

Low-rank approximation to the solution of a nonsymmetric algebraic Riccati equation from transport theory

We consider the solution of the large-scale nonsymmetric algebraic Riccati equation XCX À XD À AX þ B ¼ 0 from transport theory (Juang 1995), with M ½D; ÀC; ÀB; AŠ 2 R 2nÂ2n being a nonsingular M-matrix. In addition, A; D are rank-1 updates of diagonal matrices , with the products A À1 u; A À> u; D À1 v and D À> v computable in OðnÞ complexity, for some vectors u and v, and B, C are rank 1. The...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

ژورنال

عنوان ژورنال: Journal of Computational Physics

سال: 2023

ISSN: ['1090-2716', '0021-9991']

DOI: https://doi.org/10.1016/j.jcp.2022.111748